Cosmology with Gravitational Waves: statistical inference of the Hubble constant 
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The precise measurement of the Hubble constant Hq is among the foundations of the current 
cosmological paradigm. Due to correlations between Hq and the remaining cosmological parameters, 
a precise measurement of Hq is critical in view of future high redshift surveys. By 2015, second 
generation ground-based laser interferometers are expected to deliver a wealth of gravitational waves 
(GW) events from coalescing compact binaries up to a redshift of about 0.3. Being free of the 
systematics affecting electromagnetic measurements, GW offer the possibility of an independent 
measurement of Hq with great accuracy. This work presents a general method based on Bayesian 
inference aimed at estimating the value of the cosmological parameters for any GW event. In 
contrast to earlier work, this framework does not require the precise identification of the putative 
optical counterpart, but it considers all the potential galaxy hosts consistent with the measured sky 
position and distance posterior distributions. When applied to the planned worldwide network of 
second generation interferometers, 50 GW events will yield a measurement of Hq with an uncertainty 
of few percent. 

PACS numbers: 95.85.Sz,98.80.-k,04.30.-w 



I. INTRODUCTION 



The determination of the Hubble parameter is among 
the most critical issues in modern cosmology, with reper- 
cussions on a wide range of fields in physics ranging from 
the determination of the abundances of light elements 
in the early Universe, the content of weakly interacting 
particles and therefore about the nature of dark energy. 
While a great deal has been learnt from the original dis- 
covery of Hubble [l3|, the exact value of the expansion 
rate today is still a matter of debate. Its current ac- 
cepted value lies in the range 60 — 75 km s~^ Mpc~^[see 
[l5l , for a review] with different methods and assumptions 
yielding different values. For example, Cepheid variable 
stars in conjunction with SNe la observations yield i!fn = 
73 ± 4(statistical) ± 5(systematic) km s"^ Mpc~^ 
while the 15 years long Hubble Key Project yields 
Hq = 62.3 ± 1.3(statistical) ± 5(systematic) km s~^ 
Mpc~^. Observations of the Cosmic Microwave Back- 
ground from the Wilkinson Anisotropy Probe (WMAP) 
combined with the assumptions of a spatially flat Uni- 
verse and a constant dark energy equation of state yields 
i7o = 73 ± 3 km s"^ Mpc-^^. Without these as- 
sumptions, the data offer no additional constraints on 
Ho due to the correlation with other cosmological pa- 
rameters. The most precise measurement to date of the 
(present day) Hubble expansion rate is obtained by com- 
bining WMAP data with measurements of the Baryonic 
Acoustic Oscillations (BAO): Hq = 70.2 ± 1.4 km s"^ 
Mpc~^ [13 • Nevertheless, once one relaxes the assump- 
tion that dark energy equation of state is constant in 
time, Hq is found to lie within the range 61 and 84 km 
s~^ Mpc~^at the 2a confidence level |33|. The strong 
correlations between Hq and the other cosmological pa- 
rameters implies that uncertainties on Hq also affect es- 
timates of the other parameters. All these methods yield 



fairly narrow statistical errors but are still affected by sig- 
nificant systematic uncertainties. Observations based on 
these established techniques are expected to significantly 
clarify the picture over the next few years, but a precise 
measurement of the Hubble parameter currently remains 
an open issue, and as such blurs our understanding of the 
high redshift universe, and its mass-energy content. 

The observation of gravitational waves (GW) from co- 
alescing compact binaries potentially offers a one-step- 
only, totally independent measurement of the Hubble 
(and other cosmological) parameters, as pointed out by 
Schutz [30] over 25 years ago. In GW observations the 
luminosity distance is a direct observable [H, [26|, [30| , and 
if one could infer from other means the redshift of the 
source, one could trivially estimate the cosmological pa- 
rameters from the distance-redshift relation. As second- 
generation (or advanced) ground-based gravitational- 
wave laser interferometers are being installed, this be- 
comes a very concrete scenario, which may play a de- 
cisive role in the debate. The challenge in GW obser- 
vations is to obtain a redshift measurement for a GW 
detected source. As the GW error-box is ^ 1 — 100 deg^, 
direct redshift measurements may be very challenging, 
)ite optimistic assum] 
E, M M S and 
short-lived Gamma Ray Bursts (GRB) are associated to 
the merger of compact objects with at least a neutron 
star cornponent, this would provide such a measurement 
[m, [13, 'ssl- However, regardless of the still open debate 
on whether short GRB progenitors are indeed compact 
binaries, the fraction of coalescing systems producing a 
short GRBs might be as low as 10~^ [3]. 

These considerations bear the question of whether it is 
at all feasible to use GW binaries as a new class of stan- 
dard candles if the electro-magnetic counterpart it is not 
known. Recent literature focused on this possibility. Tay- 



despite optimistic assumptions made by several authors 
[e.g. [m, [U, [23, and many others]. For example, if 
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lor et al. [32| explored the case in which the mass function 
of neutron stars is known, while Messenger & Read [20| 
suggest a direct measurement of the redshift from GWs, 
but the required sensitivity is achievable only by third 
generation instruments. The case for space-based anten- 
nas has been investigated in MacLeod & Hogan [l9| for 
the case of Hq and in Petiteau et al. [22| for the Dark 
Energy equation of state parameter. 

In this study, I present a general method that for- 
malises in a rigorous way the original idea proposed 
in Schutz [30|. The aim is to estimate the value of 
any cosmological parameter without assuming an electro- 
magnetic identification of the source or a particular form 
of the neutron star equation of state. The results will 
be specialised to the case of ground-based gravitational 
waves detectors and I will show how they could bring 
a measurement of the Hubble parameter with good ac- 
curacy. The method is based on correlating a galaxy 
catalogue of potential hosts of the detected GW source 
with the error box in which one out of many is the pos- 
sible host and on the combination of the results from 
multiple GW detections. In fact, for second-generation 
instruments - Advanced LIGO [l6|. Advanced Virgo [5| 
and the Large Cryogenic Gravitational-wave Telescope 
[l8| are expected to provide several detections - the de- 
tection rate is estimated to be in the range ~ 1 — 100 
yr~^, depending on the actual astrophysical event rate, 
instrument duty cycles and sensitivity evolution [1] . The 
rest of the paper is organised as follows: in Section [IT] 
I describe the inference model and I specialise it to the 
case of compact binary coalescences observed from sec- 
ond generation interferometers. Section Hill presents the 
GW events catalogue on which the simulation in Section 
IIVI is based. Finally, the results are summarised and dis- 
cussed in Section IVl 



II. METHOD 

Consider a catalogue of gravitational wave events E = 
ei, . . . , e^. The posterior probability distribution for the 
cosmological parameters - the Hubble constant i^o, the 
density parameters Qrn^ ^/c, and I^a, and so on - that I 
collectively represent with the ^1, for any given event 
and cosmological model is given by: 



= J dOdzp{n,o,z\ei,n^,i) 



(1) 



where the is the set of parameters on which the GW 
waveform depends and X represents all the background 
information relevant for the problem under scrutiny. Fur- 
thermore, the integral in ([!]) is to be performed over 
the whole parameter space. For non-spinning waveforms 
= {Ai^r]^(j)c^tc^cy^S^tlj^ Dl): the chirp mass A^, the 
symmetric mass ratio 7^, the phase (j)c and time of coa- 
lescence tc^ right ascension a, declination (5, polarisation 



angle ^, inclination angle t and luminosity distance Dl- 
Thanks to the product rule, the integrand in Eq. ([1]) can 
be factorised as: 

p{n,o,z\ei,n^,i) = 

p{n\o, z, e„ n^,i)p{z\o, e,,n^,x)p(e\e,, n^,x) (2) 

where each of the factors has a simple interpretation: 
(i) p{0\ei^X)^ the joint probability distribution for the 
source parameters as inferred with parameter estima- 
tion codes [e.g. [s^, (ii) p{z\0^1-L^^X)^ the probability 
distribution for the redshift z given the observed pa- 
rameters and the assumed cosmological model and 
(iii) z, e^, the distribution for given the 

observed parameters of the source ^, the corresponding 
redshift distribution and The inference problem is 
essentially solved once one specifies appropriate distribu- 
tions. For the purpose of this paper, I will choose par- 
ticularly simple ones to exemplify the workings of this 
method. 

The posterior distribution p{0\ei.,X) is, in general, a 
complicated function of all the GW parameters and their 
correlations. Nevertheless, if one ignores these correla- 
tions it can be factorised in the product of the proba- 
bility distributions for each single parameter. Since the 
redshift distribution will depend exclusively on I)/,, a and 
(5, all the distributions for the remaining parameters will 
integrate to one. Therefore, all is left is p{Dl^ a^S\ei^X). 
I choose for this distribution a particularly simple form 
that ignores the correlations between these parameters: 



p{DL,a,S\ei,X) ex exp 



where g{a^ S) is given by: 



{Do - Dl? 
2al 



g{a,S) (3) 



g{a,5) oc 



{I 



if ^{s - 5oy 

otherwise 



cos2(<5o)(a-ao)2 < An 



(4) 

where the subscript indicates the mean values of the pa- 
rameters of interest and AQ is the uncertainty over the 
sohd angle. The value of A^} and can be heuristically 
assigned by specifying the value of the signal-to- noise ra- 
tio S/N as [i,!!^: 



Dl 
AQ 



3 X (S/N)- 



a X deg' 



f— V 



(5) 
(6) 



where a is a factor in ~ 10 — 100 that I take to be 100. 
Since I neglected all the correlations among parameters 
and these do affect the mean value I^o of the posterior for 
Dl [e.g. H^, I introduce a shift in I^o of ±/7(— 1, Ijcr^^ 
with [/(— 1,1) being a random number uniformly dis- 
tributed in the interval [—1, 1]. 

Gravitational waves alone are not able to provide 
redshift information, therefore this piece of information 
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needs to be acquired by other means [39|. Hence, past 
studies assumed a coincident GW and electromagnetic 
detection so to infer simultaneously distance and redshift 
of the source. In what follows I will relax this assump- 
tion. In this case, one can use a galaxy catalog to infer the 
possible values of the redshift associated with any given 
source. The factor ^, e^, in fact, selects 

the regions of the z space that are consistent with the 
measured 3-dimensional error box and the cosmological 
model H^. Thus, I assign: 

/ h ^kLo ^(^ - ^k) if Zk e [Zmin, ^max] and ^(tt/e, 4) < 

\ otherwise 

where Zk^o^ki ^k are the redshift and celestial coordinates 
of the /c-th galaxy and M is the total number of galaxies 
in the 3-dimensional sky box. The limits Zmin and Zmax 
are determined both by the width of the distribution for 
Dl and by the priors for ^1. For each value of allowed 
by their prior distributions, the Dl-z relation predicted 
by H^, once evaluated at the limits of the distance distri- 
bution, gives putative values of ^^min and Zmax- One can 
then choose Zmin and ^^max as the ones giving the largest 
interval. 

The probability distribution p{Q.\Dl^<^^5^ z^ei^l-L^^X) 
is instead a deterministic function. In fact, a cosmologi- 
cal model implies that the quantities Q^^Dl and z are 
not independent. For every value of the triplet I^o and 
Zk there is one function / such that /(^l, Dq^ z^) = 0. In 
particular, I^o, ^/c) = — D{Q.^ Zk) = 0. The func- 
tion D((l^z) is the luminosity distance expression as a 
function of the redshift z and of the cosmological pa- 
rameters as dictated by Hence, in full generality, 

p{^\Dl, a, S, z, e„ = S{Do - D{n, z)). 

With these choices for the relevant probability distri- 
butions, for each GW event e^, Eq. ^ becomes: 





(Do-D(a,Zfc))"l 




! d^Y.k=i! exp 


r (Do-D(Q.,Zk)Y' 
V 2aJ, J 

(8) 



with the condition that Zmin ^ Zk < ^max and g{ak^ Sk) < 
ACl. This formalism is readily extended to the full cata- 
logue of events E: 



1-n . 



vanced interferometers. The galaxy catalogue used as 
baseline for the generation of the events is the Sloan Dig- 
ital Sky Survey Data Release 7 (SDSS) [2]. The GW 
events catalogue consists of a set of 1,000 compact binary 
coalescence events distributed uniformly on the celestial 
sphere coverage of SDSS. Each galaxy is assigned equal 
probability of being the host of a GW event, therefore the 
redshift distribution of the GW follows closely the galaxy 
redshift distribution. The luminosity distance I^l is then 
calculated assuming a ACDM universe with = 70 km 
s~^ Mpc~^, Vtrn = 0-3 and I^a = 0.7. The component 
masses of the binary system are chosen by drawing from 
Ailuniform distribution mi, m2 G (1.0, IS.OjM©. As SDSS 
covers about half of the northern sky, events that are 
close to the edge of the survey coverage are not consid- 
ered. Sources were also restricted to have z < 0.1 (~ 460 
Mpc in a ACDM cosmology). This redshift corresponds 
to the putative completeness limit of SDSS and, given the 
most plausible rates [1], should yield ~ 50 detections per 
year. The full event catalogue is then divided into 20 cat- 
alogues of 50 events each. Results will be presented in the 
form of an average over these 20 independent catalogues. 
For the calculation of the S/N ^ and therefore of the un- 
certainties over distance and sky position as described by 
Eq. (O, I assumed a network of three advanced interfer- 
ometers; the two Advanced LIGOs in Hanford, WA, and 
Livingston, LA, and Advanced Virgo in Cascina, Italy. 
The sensitivity in each of the detectors is given by the 
High Power Zero Detuning configuration [28] (see Ap- 
pendix [A]). The optimal S/N ratio in each instrument is 
given by: 



S\'^ /"/isco 



df 



\h{e,f)\' 
s{f) 



(10) 



where S{f) is the power spectral density and h{0^ f) is 
the Fourier transform of the gravitational wave template. 
For h{6^f)^ I use frequency domain templates belonging 
to the so-called TaylorF2 family 0: 



HOJ)=Af 



-7/6^i^iej) 



where 



A (X 



Dl 



(11) 



(12) 



where p{Q\T-L^,X) is the prior probability distribution for 



is the amplitude of the signal from a source at lumi- 
nosity distance Dl^ and the phase ^{O^f) at the post^- 
Newtonian order is given by [J]: 



III. GRAVITATIONAL WAVE EVENTS 
CATALOGUE 

I applied the formalism presented in the previous sec- 
tion to a mock catalogue of GW events observed by ad- 
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(13) 
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In the previous expressions / is the GW frequency, 



M 
V 
M 



mi + m2 
mim2 



(14) 
(15) 
(16) 



are the total mass, the symmetric mass ratio and the 
"chirp" mass, respectively, of a binary of components 
masses rui and m2- Lastly, <l>c is the constant phase 
at some (arbitrary) reference time tc- The limits of inte- 
gration in Eq. (fTQ|) are /min = lOHz and the frequency of 
the innermost stable circular orbit: 



/isco 



(6^/V(mi+m2)) \ (17) 



The network S/N is then the rms of each detector S/N . 
Consistently with current detection searches, I considered 
only events yielding a 6'/A/' of 8 in each single interferom- 
eter. This results in a cut-off of 13.3 for the network S/N. 
Fig. [T] shows distributions for the properties of the GW 
catalogue. The network S/N ranges between 13.3 and ~ 
60, with most of the events close to threshold. This im- 
plies that typical values for the uncertainties are 25% 
on Dl and ~ 30 deg^ on the sky position. These num- 
bers are comparable with the uncertainties from second 
order Fisher Matrix estimates [35|. The typical num- 
ber of galaxies identified as potential hosts for each GW 
event ranges from a few to ^ 10,000 depending on the 
S/N ratio of the source. Note that each galaxy has the 
same probability of hosting a GW event, so the spacial 
distribution of sources follows closely the spacial distri- 
bution of galaxies. Galaxies are usually found in bound 
structures where the typical galaxy number density is 
hundred times higher than the field average, therefore 
sources are more likely to be observed from large scale 
structures such as clusters or groups of galaxies having 
approximately the same redshift. For each event, the 
most probable value for ^1 will be given by the most rep- 
resented redshift in the corresponding galaxy catalogue. 
The remaining redshifts will contribute as noise at the 
tails of the target posterior distribution. The clustering 
of galaxies was in fact the basis of the study by Schutz 
[30] . As in the aforementioned study, the combined PDF 
in Eq. (|9]) will effectively average out the tails of the pos- 
terior distribution for l1, thus allowing a more accurate 
inference over these parameters. 



IV. RESULTS 

The cosmological hypothesis considered for the 
analysis is a Friedmann- Robertson- Walker-Lemaitre uni- 
verse depending on four parameters, = h^Qrn^^k^^A 
with the condition + + = 1. In practice, the 
parameters of interest are h^^lrn^^A- As customary, I 
defined h = Hq/IOO km s~^ Mpc~^. The relevant ex- 
pression for the luminosity distance D{Cl^ z) can be found 



in prj. I choose prior probability distributions that are 
uniform in each parameter: 



h e [0.4,1.0] 

e [0.0,1.0] 
e [0.0,1.0]. 



(18) 
(19) 
(20) 



For each event, I computed the posterior distribution for 
i1, Eq. dHj), using a Nested Sampling algorithm [29[. This 
algorithm evaluates directly the denominator in Eq. (|8j), 
the evidence^ and, as a by-product, the posterior distri- 
butions of interest. I found this technique computation- 
ally more feasible that direct evaluation on a grid and, 
furthermore, more easily extendable to more general or 
complicated cosmological models. Figs. [2] and [3] show 
posterior distributions for Q for two representative events 
with low S/N and high S/N. In both cases, neither 
or l^A is measurable, their posterior distributions are es- 
sentially fiat just like their priors. For the redshift range 
under consideration {z < 0.1) the luminosity distance for 
different values of either or ft\ varies hy ^ 1%, which 
is much smaller than the typical 30% uncertainty on Dl 
given the S/N ratios. On the contrary, in both cases h is 
measurable, even though with very different accuracies. 
In the low S/N case. Fig. [21 the posterior for h tends to 
prefer values < 0.7, but with a significant support at the 
edges of the allowed range. Note that, since most of the 
potential hosts have redshifts towards the upper bound- 
ary of the allowed range from the priors on H^, the prob- 
ability distribution for h is skewed towards small values. 
A similar effect can be noted in the high S/N scenario. 
Fig. [3l even though since the 3-dimensional error box is 
very small, the posterior distribution for h is fairly well 
peaked around the correct value and it has very little sup- 
port close to the edges of the prior range. Nevertheless, 
if one would estimate h from this single event, which is 
indeed extremely optimistic, one would infer h = 0.7 ±0.1 
at the 1(7 level. The combination of the information com- 
ing from multiple sources as in Eq. ([9]), can significantly 
improve our inference on h. Fig. |4] shows the evolution 
of the joint posterior distribution for a particular reali- 
sation of 50 GW events. The comparison between the 
single event posterior, in dark gray, and the joint poste- 
rior, in light gray, highlights the effect of averaging over 
the galaxy noise at the tails of the distribution. 

h can be well determined with just few tens of obser- 
vations; ~ 25 observations will yield an estimate of the 



Hubble parameter as h = 0.70 



+0.04 
-0.05 



at 95% confidence. 



Combining 50 observations the estimate is improved to 
h 



0.71 



+0.035 
-0.035 



at 95% confidence (c/r. Fig. [5]). The 
results presented here qualitatively agree with the pre- 
diction by Schutz [30] that h is determined at the few 
percent level after 10 — 12 observations. Nevertheless, 
the uncertainty on the estimate of h does not improve 
as the square root of the number of observations, as sug- 
gested by the aforementioned study, but it is limited by 
the smallest crjjj^ , and therefore by the largest S/N ratio, 
that the instruments can achieve. 
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FIG. 1: Properties of the GW events catalogue. Top left: distribution of the relative error on Dl- Top right: distribution of 
the solid angle accuracy. Center left: distribution of recovered luminosity distances. Center right: distribution of the redshift 
of the sources. Bottom left: distribution of the network S/N . Bottom right: distribution of the number of galaxies associated 
as potential hosts to each GW event. 



V. CONCLUSIONS 

Second generation interferometers are scheduled to 
start operation in 2014-2015. Within few years of ob- 
servations, they are expected to deliver a wealth of GW 
events. GW as standard sirens offer the very realistic 
possibility of measuring Hq in way that is independent of 
the " cosmic scale ladder" traditionally used in EM-based 
cosmology. This will be a milestone in the path of un- 
derstanding the ACDM cosmological paradigm and our 
Universe in general. This paper presented a fully general 



method aimed at estimating the value of the cosmological 
parameters for any cosmological hypothesis H^, given a 
3-dimensional joint probability distribution in luminosity 
distance Dl and sky position and a probability distribu- 
tion for the redshift consistent with both the measured 
Dl and 

The application to second generation GW experiments 
shows that they can deliver a measurement of the Hubble 
parameter which is comparable with the current value de- 
rived from WMAP 7- year observations. Differently from 
most literature, this inference model makes use of all GW 
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FIG. 2: Posterior distributions for a sample GW event. The event S/N is 15.7 giving Dl = (196 ± 80) Mpc, AQ = 40 deg^ 
and a total of 1,749 potential host galaxies. The source redshift is 0.0477. Top left: distribution of the redshift of the potential 
hosts. Top right: posteriors for Qrn- Bottom left: posteriors for h. Bottom right: posteriors for Qa. 



events, without a need for the simultaneous observation 
of an EM counterpart. A scenario where one considers 
exclusively events with a putative EM detection neglects 
a significant part of the cogent information for the in- 
ference on ^1. Nevertheless, most of the information is 
gained from sources whose S/N ratio is high enough to 
provide a reasonably small error box in the sky and thus 
a peaked posterior distribution for Cl. Therefore, just one 
sure identification of the host galaxy redshift obtained, 
for example, by the simultaneous observation of a ORB 
and a GW will bring enough information to essentially 
drive our knowledge of Q. 

When compared with current EM measurements, GWs 
will not significantly improve the accuracy with which Hq 
is known. The WMAP+BAO measurement yields an un- 
certainty of ~ 2% at 68% confidence [17]. Nonetheless, 
GW and EM methods are affected by entirely different 
systematics. The simulation presented in this study high- 
lights GW as a mean of independently testing the current 
cosmological paradigm. Therefore, GW-based methods 
are to be considered not competitive, but complemen- 
tary to EM-based ones. GW are not free from systemat- 
ics. The measurement of p{Dl^ a, (5|ei,X), for example, is 
known to be affected by calibration errors that, neverthe- 
less, seem to be in general small |36|. Also the waveform 



family adopted for the analysis might introduce system- 
atic effects, although, for the low mass range considered 
in this study, waveforms are known with sufficient accu- 
racy to yield adequate measurements [7]. 

Being a proof-of-principle, this study is characterized 
by some simplistic prescriptions that, in a realistic case, 
need to be dealt with carefully. In general, the measure- 
ment of Dl cannot be disentangled from the inclination 
angle or from the sky position. This raises the issue of 
the number of detectors in operation worldwide. For in- 
stance, with only two detectors it will not be possible 
to measure Dl and sky position and thus to infer 
At the same time, if the number of detectors is greater 
than three, including LCGT [l8[ or even LIGO south 
0, the accuracy for p{Dl^ S\ei,X) is greatly improved 
[e.g. [3l|, [m and also the degeneracy between Dl and l 
is partially broken. 

The incompleteness of the survey used to obtain the 
redshifts will also affect the PDF for Q. If the "true" 
host of the GW signal is missed, the PDF for Q yields a 
low probability for the correct values, thus contributing 
essentially as noise in the the joint PDF. However, unless 
most of the single event PDFs are biased, the joint PDF 
for Q will still converge towards their true values, with 
varying convergence rate. These kind of systematics are 
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FIG. 3: Posterior distributions for a sample GW event. The event S/N is 53 giving Dl = (102 ±9) Mpc, An = 3.5 deg^ and a 
total of 120 potential host galaxies. The source redshift is 0.0218. Top left: distribution of the redshift of the potential hosts. 
Top right: posteriors for Qrn- Bottom left: posteriors for h. Bottom right: posteriors for Qa. 



especially important in preparation of the so-called third 
generation interferometers [23|. These instruments will 
observe sources in the local Universe, z < 0.1, with a 
S/N ratio ten times higher than advanced interferome- 
ters. Thus, they have the potential of providing an ex- 
tremely accurate measurement of the PDF for i^o, pos- 
sibly improving the estimate from advanced instruments 
by one order of magnitude. Furthermore, third gener- 
ation interferometers will also observe the high redshift 
Universe < 2 for neutron star binaries) providing ac- 
cess to the measurement of the energy density parame- 
ters Qrn^ and w. The analysis would further be com- 
plicated by the inclusion of the effect of weak lensing 
jiol , [l2j . Minimizing all the possible sources of bias in 
the inference of Q calls for extremely deep wide field sky 
surveys such as the planned Large Synoptic Survey Tele- 
scope (LSST) 

The framework presented in this study is fully general 
and can be applied to any GW event, regardless of the 
nature of the redshift information necessary for the infer- 
ence based on the distance-redshift relation. This might 
be provided either by deep sky surveys, as has been the 



case for this study, or by the measurement of the tidal 
couphng contribution to the GW waveform [20] or by a 
prescription based on the neutron star mass function [32| . 
The inference of ^1 using gravitational waves alone, so in- 
dependently of the "distance scale ladder", offers the very 
concrete possibility of independently testing the current 
cosmological paradigm. 
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FIG. 4: Joint posterior distribution for h after combining the single event posteriors from 1, 5, 10, 15, 20, 25, 40 and 50 events, 
respectively, for a particular realization of 50 GW events. In dark gray the single event posterior distribution corresponding 
to the latest GW event included in the calculation of the joint posterior for /i, in light gray. With increasing number of GW 
events the joint posterior becomes narrower as the noise coming from galaxies that contribute at the edges of the allowed range 
for h cancels out. 
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FIG. 5: Confidence interval evolution for Hq as a function 
of the number of events considered in the analysis. The dots 
correspond to the posterior median value obtained from 20 
realisations of 50 GW sources. The error bars correspond to 
the mean 95% confidence interval. 
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Appendix A: Advanced LIGO Zero Detuning High 
Power noise curve 



The noise curve used to calculate the S/N ratio is the 
Zero Detuning High Power design sensitivity. The ana- 
lytic fit is given by: 
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FIG. 6: The Advanced LIGO Zero Detuning High Power noise 
curve analytic fit, Eq. ()A1|) . 



